Constraining the symmetry energy at subsaturation densities using isotope binding 

energy difference and neutron skin thickness 



o 

<N 



> 
<N 

<N 
O 

m 



Zhen Zhang^ and Lie- Wen Chen*^'^ 

^ INPAC, Department of Physics and Shanghai Key Laboratory for Particle Physics and Cosmology, 
Shanghai Jiao Tong University, Shanghai 200240, China 
^Center of Theoretical Nuclear Physics, National Laboratory of Heavy Ion Accelerator, Lanzhou 730000, China 

(Dated: February 22, 2013) 

We show that the neutron skin thickness Arnp of heavy nuclei is uniquely fixed by the symmetry 
energy density slope L{p) at a subsaturation cross density pc ~ 0.11 fm~^ rather than at saturation 
density po, while the binding energy difference between a heavy isotope pair is essentially 
determined by the magnitude of the symmetry energy ^sym(p) at the same pc- Furthermore, we 
find a value of L{pc) leads to a negative EsyYn{po)-L{po) correlation while a value of ^sym(pc) leads 
to a positive one. Using data on Arnp of Sn isotopes and AE of a number of heavy isotope pairs, we 
obtain simultaneously Esym{pc) = 26.65 ±0.20 MeV and L{pc) = 46.0 ±4.5 MeV at 95% confidence 
level, whose extrapolation gives ^sym(po) = 32.3 ± 1.0 MeV and L{po) = 45.2 ± 10.0 MeV. 

PACS numbers: 21.65.Ef, 21.10.Gv, 21.30.Fe, 26.60.Gj 



The determination of density dependence of the sym- 
metry energy £^sym(p), which characterizes the isospin 
dependent part of the equation of state (EOS) of asym- 
metric nuclear matter, is of fundamental importance due 
to its multifaceted roles in nuclear physics and astro- 
physics [1-4] as well as some issues of new physics be- 
yond the standard model [El-Q- Due to the particularity 
of nuclear saturation density po (~ 0.16 fm~^), a lot of 
works have been devoted to constraining quantitatively 
the magnitude and density slope of the symmetry energy 
at Po, i-e., Esym{po) and L{po), by analyzing terrestrial 
nuclear experiments and astrophysical observations. Al- 
though significant progress has been made during the last 
decade, large uncertainties on ^sym(po) and L{po) still 
exist (See, e.g., Refs. [ifl, [Mil). For example, while 
the value of £^sym(po) is determined to be around 30 ± 4 
MeV, the extracted L{po) varies drastically from about 
20 to 115 MeV, depending on the observables and analy- 
sis methods. To better understand the model dependence 
of the constraints and reduce the uncertainties is thus of 
critical importance and remains a big challenge. 

Among the experimental observables, the neutron skin 
thickness Ar^p = (^n)"^^^ ~ (^p)"^^^ of heavy nuclei, i.e., 
the difference of the neutron and proton rms radii, has 
been shown to be a good probe of £^sym(p) [l2l-[22|, and 
this provides a strong motivation for the Lead Radius Ex- 
periment (PREX) being performed at the Jefferson Lab- 
oratory to determine the (r^)"^/^ of ^°^Pb to about 1% 
accuracy by measuring the parity-violating electroweak 
asymmetry i n the e lastic scattering of polarized electrons 
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sion [l3| . Since the pressure of neutron-rich matter is es- 
sentially controlled by the density dependence of £^sym(p) 
and the characteristic (average) density in finite nuclei is 
less than po (See, e.g., Ref. ^]), one expects that the 
Avnp should depend on the subsaturation density behav- 
iors of the Esym{p) 0. Brown and Typel [lij noted firstly 
that the Avnp of heavy nuclei from model calculations 
is linearly correlated with the pressure of pure neutron 
matter at a subsaturation density of 0.1 fm~^. The lin- 
ear correlation of the Avnp with both ^sym(po) and L 



has also been observed in mean-field calculations |12l-l2C 
using many existing nuclear effective interactions. 

Recently, a remarkable negative correlation between 
^sym(po) and L{po) has been obtained by analyzing ex- 



from ^^^Pb |23l-[25j. Physically, the Ar^p depends on 
the pressure of neutron-rich matter in nuclei which will 
balance against the pressure due to nuclear surface ten- 



isting data on Ar^p of Sn isotopes [27|, showing a striking 
contrast to other constraints that essentially give a pos- 
itive Esym{po)-L{po) correlation (See, e.g., Refs. |9l4l l\ ) . 
It is thus of great interest to understand the physics be- 
hind this negative £^sym(po)-^(po) correlation. Using a 
recently developed correlation analysis method, we show 
here that the Avnp of heavy nuclei is uniquely fixed by 
the L{p) at a subsaturation cross density pc ^ 0.11 fm~^, 
which naturally leads to a negative £^sym(po)-^(po) cor- 
relation. Furthermore, we demonstrate that the bind- 
ing energy difference between a heavy isotope pair is 
essentially determined by the magnitude of £^sym(p) at 
the same pc- The isotope binding energy difference and 
neutron skin thickness of heavy nuclei can thus be used 
to constrain simultaneously the magnitude and density 
slope of the £^sym(p) at pc ^ 0.11 fm~^. 

The EOS of asymmetric nuclear matter at baryon den- 
sity p and isospin asymmetry S = {pn — Pp) / {pp ^ pn) ^ 
given by its binding energy per nucleon, can be expanded 
to 2nd-order in 5 as 



E{p, 8) = Eo{p) ± E,y^{p)S^ ± 0((5^), 



(1) 



* Corresponding author: lwchen@sjtu.edu.cn 



where Eo{p) = E{p^6 = 0) is the binding energy per 
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nucleon in symmetric nuclear matter, and the nuclear 
symmetry energy is expressed as 



^sym(p) = 



1 d^E{p,5), 



\s=o- 



(2) 



2! 85^ 

Around a reference density p^, the £^sym(p) can be char- 
acterized by using the value of £^sym(Pr) and the density 
slope parameter L{pr) = 3pr 



a{p) I 



^sym 



■^sym 



(Pr) 



dp 
L{pr) 



Xr + 0(Xr), 



(3) 



with Xr = {p- pr)/pr' 

In the present work, we use the Skyrme-Hartree-Fock 
(SHF) approach with the so-called standard form of 
Skyrme force (see, e.g., Ref. [28]) which includes 10 pa- 
rameters, i.e., the 9 Skyrme force parameters <j, to — ^3, 
^0 — ^3, and the spin-orbit coupling constant Wo- This 
standard SHF approach has been shown to be very suc- 
cessful in describing the structure of finite nuclei, es- 
pecially glob al prop erties such as binding energies and 



charge radii [28 



-30|. Instead of using directly the 9 
Skyrme force parameters, we can express them explicitly 
in terms of 9 macroscopic quantities, i.e., po, £^o(po), the 



incompressibility Kq^ the isoscalar effective mass m 



the isovector effective mass m* o, 



s,05 

EsymiPr), L{pr), the 

gradient coefficient G^, and the symmetry-gradient coef- 
ficient Gv' Then, by varying individually these macro- 
scopic quantities within their known ranges, we can ex- 
amine more transparently the correlation of properties of 
finite nuclei with each individual macroscopic quantity. 
Recently, this correlation analysis method has been suc- 
cessfully applied to study the neutron skin [27[ and giant 
monopole resonance of finite nuclei [3l| with p^ = po. 

To examine the correlation of the Ar^p of heavy nuclei 
with each macroscopic quantity, especially on £^sym(Pr) 
and L{pr), we show in Fig.fflthe Ar^p of ^^^Pb from SHF 
with the Skyrme force MSLO ^3 t>y varying individually 

L{pr), Gv, Gs, Eq{po), ^sym(pr), Kq, mJo, ^^,0^ PO, 

and Wo within their empirical uncertain ranges, namely, 
varying one quantity at a time while keeping all others 
at their default values in MSLO, for three values of pr, 
i.e., pr = 0.06, 0.11, and 0.16 fm~^. It is seen from 
Fig.[T]that the Avnp of ^^^Pb exhibits a strong correlation 
with L{pr) and £^sym(Pr) while much weak correlation 
with other macroscopic quantities. In particular, while 
the Arnp always increases with L{pr), it can increase 
or decrease with £^sym(pr), depending on the value of 
pr. Most interestingly, one can see that, for p^ = p^ ^ 
0.11 fm~^, the Arnp becomes essentially independent of 
Esymipr) and only sensitive to the value of L{pr). These 
features imply that the L{pc) is a unique quantity to 
determine the Avnp of heavy nuclei, and the experimental 
data on Arnp of heavy nuclei can put strong limit on the 
value of L{pc). 

In order to determine the magnitude of the symmetry 
energy at i-c, Esym{Pc), we propose here to use the 
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FIG. 1: (Color online) The Arnp of ^^^Pb from SHF with 
MSLO by varying individually L{pr) (a), Gv (b), Gs (c), 
Eo{po) (d), Esymipr) (e), Ko (f), m*,o (g), ml^o (h), po (i), and 
Wo (j) for Pr = 0.06, 0.11, and 0.16 fm"^ The Esym(pr) is 
shifted by adding 15 and 7 MeV for pr = 0.06 and 0.11 fm"^, 
respectively. 
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FIG. 2: (Color online) Same as Fig. [T] but for between 
isotope pair ^^^Pb and ^'^^Pb. 



difference of binding energy per nucleon between an iso- 
tope pair, denoted as AE". From the well-known semiem- 
pirical nuclear mass formula, one can easily see that the 
AE essentially reflects the symmetry energy of finite nu- 
clei with mass number A, i.e., cisyin 

(A), especially for the 
spherical even-even nuclei. For heavy nuclei , on e has the 
empirical relation of asym(^) ^ ^sym(Pc) [S, EH, and 
thus expects that the AE for heavy isotope pairs should 
be a good probe of Esym(Pc)- Similarly as in Fig. (TJ 
we show in Fig. [2] the AE for isotope pair ^^^Pb and 
^^^Pb from SHF with MSLO by varying individually the 
9 macroscopic quantities and Wo for pr = 0.06, 0.11 and 
0.16 fm~^. As expected, it is seen from Fig. [2] that, for 
pr = 0.11 fm~^, the AE exhibits a very strong corre- 
lation with Esymipr) while it displays relatively weak 
dependence on L{pr), Gv, Gs, ml q and Wq, and es- 
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FIG. 3: (Color online) ^ function of (pc) (a) and 

L(pc) (b). Dashed lines represent results from SHF with 
MSLO while the solid circles are obtained by optimization. 



sentially no dependence on other macroscopic quantities. 
For pr = 0.06 and 0.16 fm~^, the AE also exhibits strong 
correlation with L{pr). These results indicate that the 
AE for heavy isotope pairs indeed provides a good probe 
of Esym{pc)' We note here that using other standard 
Skyrme forces instead of MSLO or using other heavy nu- 
clei such as Sn isotopes leads to the similar results as 
shown in Fig. [T] and Fig. [2l We have also checked with a 
number of existing standard Skyrme forces and obtained 
the similar conclusion (33| . 

Experimentally, there are very rich and accurate data 
on ground state binding energy of finite nuclei. To 
constrain £^sym(Pc) from AE^ here we select 19 heavy 
isotope pairs which are all spherical even-even nu- 

1 . 1 218,206t) 216,194t3 214,178t3U 212,178t3U 

clei, namely, sfi^^^ 84 82^^)7 82 ^t), 

210,178p^ 208,178p^ 206,178p^ 206,17Tjjg 136,l(5Brp^ 
132,l§§ci_ ' 132,l§lo_' 132,l§io_ ' 132,l§8o_ ' 132, l5§ ' 



Sn, 



50 



Sn, 



50 



Sn, 



50 



Sn, 



5QSn, rQk 

132,114q^ 130,98^-1 124,96r);j" 94,84^;r^ "^^^^^ 94,82r7^ 

5oSn, \s^^^ 46^^^ 42^0, and \qZt. On 
the other hand, the Ar^p of heavy Sn isotopes has been 
systematically measured using various methods, and here 
we use the existing 21 data on Ar^p of Sn isotopes 33-[39| 
to constrain the L{pc). Firstly, with all other parameters 
fixed at their default values in MSLO, we calculate the 
from the difference between theoretical and experimen- 
tal AE {Arnp) with different Esym{pc) {L{pc))^ and the 
results are shown by dashed lines in Fig. [3l When we 
evaluate the we use the experimental errors for Arnp 
while assign a theoretical relative error of 23% for AE 
so that the minimum value of is close to the number 
of data points to make the x^-analysis valid [ioj. From 
the dashed lines in Fig. [3] obtained with MSLO, we can 
extract a value of £^sym(Pc) = 26.08 ± 0.17 MeV and 
= 47.3 ± 4.5 MeV at 95% confidence level. 
The values of ^sym(Pc) = 26.08 ± 0.17 MeV and 
L{pc) = 47.3 ± 4.5 MeV are obtained by assuming 
the other macroscopic quantities are fixed at their de- 



fault values in MSLO when £^sym(Pc) or L{pc) is var- 
ied, and thus the correlations of AE {Arnp) with the 
other macroscopic quantities are totally neglected. To 
be more precise, one should consider the possible vari- 
ations of the other macroscopic quantities due to the 
correlations introduced by fitting the theoretical predic- 
tions to some well-known experimental observables or 
empirical values. Taking this into account, for a fixed 
Esymipc) {L{pc))^ we optimize all the other parameters, 
instead of simply keeping them at their default values 
in MSLO, by minimizing the weighted sum of evalu- 
ated from the difference between SHF prediction and the 
experimental data for some selected observables using 
the simulated annealing method [4l|. In the optimiza- 
tion, we select the binding energy per nucleon Eb and 
charge rms radii Tc of 25 spherical even-even nuclei, i.e., 

204, 202, 200, 198, 196, 194, 192, 190pT 124,122,120,118,116,112,108 q 

82 ' 50 ' 

64,62,60,58^T- 50,48,46,44,42,40^ a ,1 n 

'23iNi, and 20^^- ^ theoretical error 

of 0.013 MeV is assigned to Eb and 0.011 fm to Tc so 
that the respective is roughly equal to the number 
of data points. We further constrain the macroscopic 
parameters in the optimization by requiring (1) the neu- 
tron 3pi/2 — 3p3/2 energy level splitting in ^^^Pb should 
lie in the range of 0.8 — 1.0 MeV; (2) the pressure of sym- 
metric nuclear matter and the binding energy of pure 
neutron matter should be consistent with constraints ob- 
tained from flow data in heavy ion collisions [42| and 



the latest chiral effective field theory calculations with 
controlled uncertainties [isj, respectively; (3) the critical 
density pcr, which is determined by the stability condi- 
tions from Landau parameters, should be greater than 
2po; and (4) rrig Q should be greater than m* q and here 
we set m* Q — ml q = 0.1m (m is nucleon mass in vacuum) 
to be consistent with the extraction from global nucleon 
optical potentials constrained by world data on nucleon- 
nucleus and (p,n) charge-exchange reactions [44]. In the 
optimization, for a fixed £^sym(Pc) (^(Pc)), the contribu- 
tion of the 21 data on Arnp of Sn isotopes (the 19 AE" 
data) is included in the weighted sum of x^- 

By using the optimized values for other macroscopic 
quantities at a fixed £^sym(Pc) (^(Pc)), we show the x^ 
evaluated from the difference between theoretical and ex- 
perimental AE (Arnp) as a function of £^sym(Pc) {L{pc)) 
by solid circles in Fig. [H From the new relation be- 
tween x^ and Esymipc) {L{pc))^ we can extract a value of 
Esymipc) = 26.65±0.20 MeV and L{pc) = 46.0±4.5 MeV 
at 95% confidence level. It is seen that the Esym{pc) = 
26.65 ± 0.20 MeV obtained with optimization is larger 
than the £^sym(Pc) = 26.08 ± 0.17 MeV extracted using 
MSLO without optimization by about 0.6 MeV, imply- 
ing the correlations of AE" for heavy isotope pairs with 
L{pr)^ Gy, G^, mJo7 and Wq observed in Fig. [2] play a 
certain roles in the extraction of Esym(Pc)- On the other 
hand, it is remarkable to see that the extracted value of 
L{pc) = 46.0 ± 4.5 MeV with optimization agrees well 
with L{pc) = 47.3 ±4.5 MeV obtained using MSLO with- 
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TABLE I: Skyrme parameters in MSLl (left side) and some 
corresponding nuclear properties (right side). 
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out optimization, indicating the correlations of Avnp of 
heavy nuclei with other macroscopic quantities play mi- 
nor roles in the extraction of L{pc), as shown in Fig. [TJ 
Furthermore, by including simultaneously the 21 data on 
Arnp of Sn isotopes and the 19 data of AE in the opti- 
mization, we obtain a globally optimized parameter set 
of Skyrme force denoted as MSLl as shown in Table HI 
As expected, the MSLl nicely reproduces the optimized 
values of £^sym(Pc) and L{pc) obtained from analy- 
sis with optimization for AE and Ar^p, respectively, as 
shown by solid circles in Fig. [3l 

The extrapolation is necessary to obtain information 
on the symmetry energy at po from £^sym(Pc) = 26.65 ± 
0.20 MeV and L{pc) = 46.0 ± 4.5 MeV. Shown in Fig. 1 
are contours in the Esym{po)-L{po) plane for £^sym(Pc) 
and L{pc) from SHF calculation with MSLl. It is in- 
teresting to see that a value of L{pc) leads to a neg- 
ative Esym{po)-L{po) correlation, which was observed 
previously by directly constraining ^sym(po) and L{po) 
from analyzing the Ar^p data of Sn isotopes On 
the other hand, a value of Esym{pc) leads to a positive 
Esym{po)-L{po) correlation as expected [32|. Combing 



FIG. 4: (Color online) Contour curves in the Esym{po)-L{po) 
plane for Esym{pc) and L{pc) from SHF calculation with 
MSLl. The shaded region represents the overlap of con- 
straints obtained from -Esym (pc) = 26.65 ± 0.20 MeV and 
L{pc) = 46.0 ± 4.5 MeV. 



difference AE between a heavy isotope pair is essentially 
determined by the magnitude of the symmetry energy 
^sym(p) at a subsaturation cross density pc ~ 0.11 fm~^, 
the symmetry energy density slope L{p) at the same pc 
fixes uniquely the neutron skin thickness Avnp of heavy 
nuclei. Furthermore, we find a value of L{pc) leads to 
an anti-correlation between £^sym(po) and L{po) at sat- 
uration density po while a value of £^sym(pc) leads to a 
positive Esym{po)-L{po) correlation. The existing data 
on Arnp of Sn isotopes and AE for a number of heavy 
isotope pairs put simultaneously stringent constraints on 
the magnitude and density slope of the £^sym(p) at pc^ 
i.e., Esymipc) = 26.65 ±0.20 MeV and L{pc) = 46.0 ±4.5 
MeV at 95% confidence level, whose extrapolation gives 
the constraints from E^ymipc) = 26.65 ± 0.20 MeV and £^sym(po) = 32.3 ± 1.0 MeV and L{po) = 45.2 ± 10.0 



L{pc) = 46.0 ± 4.5 MeV leads to a quite stringent con- 
straint on both Esym{po) and L{po) as indicated by the 
shaded region in Fig.lH giving £^sym(po) = 32. 3 ±1.0 MeV 
and L{po) = 45.2 ± 10.0 MeV, which is essentially con- 
sistent with other constraints extracted from terrestrial 
experiments and astrophysical observations as well as 
theoretical calculations with controlled uncertainties ^ 
11, [4^ but with much higher precision. Furthermore, 



we find the extracted £^sym(Pc) = 26.65 ± 0.20 MeV and 
L{pc) = 46.0 ± 4.5 MeV in this work leads to a quite 
strong constraint of Ar^p = 0.170 ± 0.016 fm for ^^^Pb 
and pt = 0.082 ± 0.005 fm~^ for the core-crust transi- 
tion density of neutron stars evaluated in a dynamical 
approach [45[. 

In summary, we show that while the binding energy 



MeV. The extracted £^sym(Pc) and L{pc) also leads to a 
strong constraint of Avnp = 0.170 ± 0.016 fm for ^^^Pb 
and pt = 0.082 ±0.005 fm~^ for the core-crust transition 
density of neutron stars. 
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